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We have studied the effects of the Standard-Model Extension (SME) on hydrogen as a realization 
of new physics effects that incorporate Lorentz and CPT violation. Specifically, we calculated the 
SME-induced energy level shifts at order a 2 times the SME parameters. We obtained contributions 
at this order both from the non-relativistic effective Hamiltonian for motion of a spin- 1/2 particle 
in the presence of SME interactions and also from SME corrections to the propagator for exchange 
photons. We applied our result to the 2S — IS transition in hydrogen, which has been measured 
with extremely high precision. The results obtained in this work give the leading SME corrections 
for this transition. 
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I. INTRODUCTION 



For the past century atomic hydrogen has played a unique role at the intersection between quantum theory and 
experiment. One need only think of the Bohr model, Schrodinger's wave mechanics, hydrogen fine structure, the 
fS | Lamb shift, the hydrogen maser, hydrogen condensates, the importance of hydrogen in astrophysics and cosmology, 
' the multitude of ultra precise spectroscopic measurements, and hydrogen's use in the least squares adjustment of 
the fundamental constants [l[ (to name a few) in order to recognize the central and enduring status of hydrogen in 
<-j '. physics. Because of the well-developed theory and precise measurements hydrogen is an attractive test bed for 

'— '« ruling out or setting bounds on new physics. In this work we will extend the study of possible violations of Lorentz 
and CPT symmetry in hydrogen in the context of the Standard-Model Extension (SME) of Colladay and Kostelecky 

The SME is a general framework for describing the effects that small violations of Lorentz and CPT symmetries 
would have on standard model processes. It is based on the idea that Lorentz and CPT violation might arise through 
spontaneous symmetry violation of a more fundamental theory [7J. The phenomenological effect of the SME is 
C*~) , to add a number of new interactions to the standard model formed from the usual standard model fields coupled 
1 to spacetime constant background tensor fields. These background fields are operationally just sets of symmetry- 
t-H I breaking parameters organized as scalars, vectors, and tensors of various ranks and symmetries. The implications of 
the SME have been explored in many areas of physics. A systematic overview organized as an extensive set of tables 
of bounds on the SME coefficients has been given by Kostelecky and Russell @ . Reviews of the SME have been given 
'. in Refs. (§£□]]. 

The SME causes small shifts in the energy levels of hydrogen that could, in principle, be detected by high precision 
spectroscopy. In standard theory, the Dirac equation with a Coulomb potential is a convenient starting point for 
obtaining the energy levels of hydrogen-giving levels correctly through fine structure at the O(rn e o: 4 ) level. The SME 
leads to a modified Dirac equation, and through a non-relativistic reduction, to a modified effective Hamiltonian for 
use in a standard quantum mechanical framework. The corrections to the effective Hamiltonian have the form of an 
expansion in powers of the relative momentum. Since expectation values of the relative momentum p have a typical 
order m e a (in units where h = c= 1) where a ~ 1/137.036 is the fine structure constant, the low order terms in this 
expansion can be used to find the leading perturbations to the energy levels. The leading SME corrections to hydrogen 
energy levels were worked out by Bluhm, Kostelecky, and Russell [l2[ and Shore and additional leading order 
effects in hydrogen were considered by Ferreira and Moucherek [l4| and Kharlanov and Zhukovsky [l5[ . By "leading 
order" we mean those corrections linear in SME parameters and of zeroth order in a. We will never consider terms 
of higher order than linear in the SME parameters since those should be completely negligible. The "higher order" 
corrections that we obtain in this work are higher order in a. Results at order a 2 times a subset of SME parameters 



have been reported by Bluhm et al. [121 ] and by Altschul [16] : our results are more general. Other Lorentz- violating 



corrections in hydrogen were analyzed by Belich et al. [Vi 



* gadkins@fandm.edu 



2 



In this work we obtain the corrections to the hydrogen energy levels at order a 2 times the SME parameters for 
SME corrections to the electron, proton, and photon interactions in hydrogen. We start with the SME-corrected 
Dirac equation and work out its non-relativistic expansion. Our result agrees with that obtained earlier [l8l - l20l |. We 
work out all corrections coming from terms linear in p (which vanish) and quadratic in p in the effective Hamiltonian. 
The new corrections are at order a 2 times SME parameters. We obtain the complete set of such corrections in a 
compact form. Next, we calculate the contributions at this order coming from SME corrections to the propagators 
of exchanged photons. Finally, we apply our result to the 25 — IS transition in hydrogen, which has been measured 
with particularly high precision [21] . 

II. EFFECTIVE HAMILTONIAN FOR LOW-ENERGY SME CORRECTIONS 

As a first step in obtaining the SME effects on non-relativistic systems, we find the extended Dirac equation for a 
fermion of mass m in the presence of an electromagnetic field. The appropriate SME Lagrangian has the form [22j 

L = h,^Y v D v ^-^M^ (1) 

where D v = d v + iqA v is the covariant derivative, — d^Ay — d v A^ and T v and M contain the SME coefficients: 

M = m + a^ + b^^ + ^H^a^. (2) 

All of the SME coefficients are real; c pi/ and g? m „ are traceless, is antisymmetric, and gx^ is antisymmetric on the 
A/i indices. Our gamma matrices are defined using the conventions of Itzykson and Zuber [23j and we use a timelike 
metric. While all of the SME parameters are Lorentz violating, the CPT-odd parameters a M , b^, e M , / M , and g x ^ v 
also govern CPT violation [(| [22j ■ The Dirac equation that results from this Lagrangian is 

(iT^D^ - M) $ = 0. (3) 

The naive Hamiltonian associated with this Dirac equation is found by rearranging ([3]) into 

iT°doiJ = (-»f • V + gTM„ + M) V, (4) 

or idgip = Hip where 

H = (r°) _1 (-if • V + gP'A, + M) (5) 



However, H is not hermitian and so is unacceptable. As shown by |18l - l20| . we can repair the hermiticity problem by 
first making a spacetime-constant spinor field redefinition in the original Lagrangian: 

^ - A X , (6) 
where A is chosen to restore the usual time derivative coupling. In particular, we choose 

A= ( 7 r°r 1/2 . (7) 

This matrix is positive definite, invertible, and hermitian. The modified Dirac equation takes the standard form 
i<9oX = Hx where the hermitian Hamiltonian is 

H = A-/°f ■ ttA + Aj°MA + qA° (8) 

with 7? = p — qA. We expand A to first order in the SME parameters and then apply the Foldy-Wouthuysen procedure 
[24| to obtain the nonrelativistic expansion of the Hamiltonian. We apply a succession of unitary transformations that 
act to diagonalize the Hamiltonian, using powers of m to keep track of the order since each 1/m is associated with a 
factor of 7? and if is small in a non-relativistic system. We specialize to the situation where ^4° is the Coulomb potential 
and A is a time- independent vector potential that leads to a constant magnetic field B. We work to first order in the 
magnetic field and neglect corrections to it involving either a or the SME parameters. We have retained terms giving 
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non-SME corrections of order ma 4 , the term of first order in B, and all terms of 0(SME ■ a°), 0(SME ■ a 1 ), and 
0(SME - a 2 ). The effective fermion Hamiltonian at this level of approximation is 

+ (A + B k a k ) + (C + D lk a k ) £ + (£y + ^ fc a fe ) + 0(SM£ • a 3 ) (9) 

where 

A = ciq — meo — tocoo, (10a) 

Bk = -bk + md k0 - ^e ka b (mg a bo - H ab ) , (10b) 

d = <n —m(c i0 + coi) —mei, (10c) 

Di k = -b 5 ik + e ika H a o + m (d ki + d oa S lk ) - -me a bkg a bi - mei ka ga00, (lOd) 

Eij = -mcij - -mcooSij, (lOe) 

Fijk = S jk di + -e ka {taipbp - 2mg a0l - mg ai0 ) , (lOf) 

with di = mdoi + ^mdiQ — \ti a bBab- This effective Hamiltonian is a sum of the usual non-relativistic Hamiltonian 

Hnr = m + p 2 /(2m) + qAo, a magnetic a ■ B interaction, fine structure contributions (relativistic kinetic energy, 
spin-orbit, Darwin), and the SME corrections. We write 0(a 5 ) and 0(SME ■ a 3 ) to indicate the orders of omitted 
terms when the fermion is bound in a Coulomb field. While we distinguish between contravariant and covariant 
indices on SME tensors a M , b^, etc., we do not distinguish between superscripted and subscripted indices on purely 
three-dimensional objects like the Pauli matrices <Tj = a 1 , <5y , eyt, and the coefficients B k , Cik, ■ ■ ■ , Fijk- We note 
that ao and eo enter only as a uniform shift to all levels, and / M drops out entirely to this order (as expected: see 
[25M27I ] ) . Our expression for the effective Hamiltonian is in agreement with the known result [l8l - [20l | . 



III. FERMION-SECTOR ENERGY SHIFTS 



In this section we will work out the energy shifts on hydrogen states due to the SME effects on the electron and 
proton up to corrections of 0(SME - a 2 ). We first describe the states in the absence of the SME. The electron 
variables are labeled by n, I, j, rrij where the total electron angular momentum J = L + S has a quantum number 
j — £ ±1/2 in terms of the orbital angular momentum quantum number £. The states |n, £,j, rrij) are approximate 
eigenstates of the Hamiltonian including fine structure and Lamb shift contributions. The small corrections to these 
states that appear in the exact eigenstates (i.e. including fine structure and Lamb shift) do not affect SME corrections 
to 0(SME ■ a 2 ). Explicit forms for the \n,l,j, rrij) states are 

\n,£,j = £±1/2,7713) = ±c{£,±m J )\n,£,m j - 1/2) | + l/2) e + c(l, ^mj)\n,t., rrij + 1/2) | - l/2) e , (11) 

where the \n,£,m) are the spinless Coulomb bound states with (f\n,£,m) = R n i(r)Y/! n (9, 4>), the \m s ) e are the 
electron spin states, and 

x f£ + m + l/2\ 1/2 , x 

More explicitly, one has 

^=W. mj (?) = (r\n, £, j = £ ± 1/2, mj > = i^(r) ( ±f ff ' ^ 3 ^J/2^) ■ (13) 

y c{£,Tmj)Y t (d,4>) J 

We abbreviate these states as \n,£,j =£± 1/2, rrij) = \mj)±. Only the + states exist when £ = 0. 

The proton spin / affects the energy levels through the magnetic hyperfine interaction. The total angular momentum 
F = J + 1 = (L + S) + 1 has quantum number / = j ±' 1/2. It is conserved due to rotational symmetry, so the states 
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having various values of F z , with eigenvalues hint, would be degenerate in the absence of an external magnetic field, 
but all other degeneracies are broken. These states have the form 

\n,l,j = £±l/2,f=j±'l/2,m J ) = ±'c(j, ±'m f )\m f - 1/2 }± | + l/2) p 

+c{j,T'm f )\m f + l/2)± \ - l/2) p 
= \m f ) ±t± , (14) 

where the \mi) p are the proton spin states. We imagine our atom to be in the presence of a small external magnetic 
field so that the states |to/)±,±' are all non-degenerate. This allows us to use the formalism of non-degenerate 
perturbation theory to calculation the SME corrections. In section fVTl we will discuss the situation where the external 
magnetic field is not present. 

Both the electron and the proton in the hydrogen atom are affected by SME corrections of the form given in © , so 
there should really be one set of SME corrections involving the electron variables (mass, momentum, spin matrices, 
and SME coefficients) and a second set of SME corrections involving the proton variables. The electron and proton 
momenta sum to zero (in the center of mass frame), and have magnitude proportional to \ia « m e a in a Coulombic 
bound state since the reduced mass fi « m e . It would seem that proton corrections such as p % p> /m^ would be 
negligible compared to electron ones from p % pP /m 2 given the small mass ratio m e /m p w 1/1836, but we have no a 
priori knowledge of the size of the electron SME coefficients relative to the proton SME coefficients, so corrections 
of order 0(SME ■ a 2 ) for the proton might be comparable in size to those for the electron. In any case, we have 
calculated both electron and proton corrections to order 0(SME ■ a 2 ). 

We first consider the 0(SME ■ a ) corrections that arise from terms in H e g that are linear in the SME coefficients 
and independent of p. These terms are spatial constants so the matrix elements are easy to evaluate. We can write 
these terms as 

H cS -> A e + AP + B e k <j k + B p k a k (15) 

with A e and A p given by expressions like (|10ap and Bf. and B v k by expressions like (jlObj) . We use the Wigner-Ekhart 
theorem [28| to relate the / = 1 operators a e and a p to F inside the expectation values, finding t? e — > £, e F, a p — > £ p F 
where the explicit values 

_ 2(2j + 1 T' 1) , 2 
* e (2^+l)(2j + l)' * p_ 2j + l [ ' 

were obtained by direct evaluation of the expectation values using the states (fl4)l . We find the 0(SME ■ a ) energy 
shift to be 

£ S °me = {A e + AP + Bla k e + Bi<j k p ) |m / )±,±, = (A e + A*) + + &i?£) m f . (17) 

This z-axis is determined by the direction of the magnetic field and will be different from the Z-axis of the conventional 
sun-centered inertial frame [1^, [3(| ■ 

The terms in H e g that are linear in the momentum make contributions of 0(SME - a 1 ) to matrix elements, but do 
not contribute to expectation values. One way to understand this is by thinking about parity: the states |m/)± i ±' of 
(|14p have definite orbital angular momentum I and parity (— 1) , and the matrix element of p (an odd-parity operator) 
between two states of the same parity must vanish. 

The terms in H c g that are quadratic in p make contributions of 0(SME ■ a 2 ). These terms have the forms 



pipP p^p^u k 

H e g -> EL — 2" + F ijk T~ ( 18 ) 

mi m 



for the electron, with corresponding terms for the proton. The energy corrections are given by the expectation value of 
(|18p in the states (fl"4"]) . so we must evaluate (m/lpV 3 'o"e l m /)±,±' (where a® = 1) in these states. We find it convenient 
to evaluate these expectation values in momentum space using wave functions 

*.>«->/Ws - ft-w (19) 

instead of p3p . where now R n i(p) is the momentum space radial function (given in terms of Gegenbauer polynomials, 
see [3l|, HH) and 8, <fi represent the angles of p. The radial dependence of these expectation values factors out: 

(m/|pV'o'e|m/}± ) ±' = (m f \p 2 p l p j a*\m f )±. ± , = (p 2 )n(m f \p i p> a^\m f )±^ . (20) 
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TABLE I: Table of coefficients for writing p l p J in terms of Y^iO, <f>). The coefficient = equals y / 2n/15 times the number 
given at the appropriate spot in the table. 
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The expectation value of p 2 can be evaluated using the virial theorem in the usual way: 

1 2 

(p 2 ) n = -2m e E n = S-— (21) 
where E n = —m e a 2 /(2n 2 ) are the Bohr energies. We decompose p % p> into irreducible tensors by writing 

ffP = X tj + ~£y (22) 

where 

i 2 

X<i = - & = E C «YfV> </>)■ ( 23 ) 

The coefficients are given in Table HI 

We build up to the evaluation of (m,f\p t pi<T£ |m/)±,±' in three steps. First, we obtain the matrix elements of p % jp 
in the spinless states \n,£,m). We use these to find the expectation values {mj\p % pio*\mj)± in the states \rrij)± of 
the electron with spin. Finally we evaluate the expectation values (mf\p % p'a% |m/)±,±/. 

The spinless matrix elements are easy to obtain by use of the decomposition of described in (122|l and ()23|) . 
One has 

(ri, I, m|py |n, I, m') = (n, m\ j YJ ^"(fl, 0) + ^ j |n, rri) 



where the integral over solid angle 



<mYr(e^)YT- m '{o^)Yp\e^) = ^cj^^c^ (25) 



is given in terms of Clebsch- Gordon coefficients [28[. Matrix elements of p l pi are obtained from those of p l pP by 
multiplication by (p 2 ) n - 

When electron spin is included we need the matrix elements 



KW<K)±= J dn (±c(£,± mj )Yp 1/2 (9^),c(£,Tm J )Yp +1/2 (e,cj ) )yfpia 



, (±c{£,±m' j )Y< 1/2 (6 
c(£,Tm'AY^ +1/2 (9, 



The 2x2 matrix algebra gives, for example, 



(26) 



(m j \p i p>a° e \m' j ) ± = c{l,±m j )c{l,±m' i ){n,l,m i - l/2\p l p>\n,£,m' 1 - 1/2) 



j / ± ' — " J - J -" l 'jy l -'\ > J jJ \ ' ' 3 I Mr 1 F | 9 5 

+ c{£,^mj)c{£,Tm'Mn,£,mj + l/2\p l p>\n,£,m' 1 + 1/2), (27) 
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with analogous expressions when k is 1, 2, or 3. It turns out that we need diagonal matrix elements only. The results 
are: 

KIpYK> ± = ^ 2i(j - + 1) J ) 2 + ( 2j -0- + l) J ^ (28a) 

... / (2£+l) 2 -4m 2 \ 

K|p*pV e K>± = T m, ^ 2(2 , + 3 )(2<+ 1)(2 j_ 1} ) (ft^a + ftafja) , (28c) 

(m-W^\m-) ± - ±m 3 )(2l - 1) T 2(21+ 1) + 4m 2 \ ^ + ^ 

^ipp^lm^i - ±m^ (2* + 3)(2* + l)(2*-l) J 2 

^ (21 + 3)(2l - 1) ±2(21+1)- 4m 2 ^ 
± ( (2l + 3)(2*+l)(2*-l) J ^ (28d) 

Finally, we are in a position to obtain the full expectation values (mf\p l p , a£ |to/)±.±<. We find these expectation 
values by using the states (|14[) and noticing that the operator p l p^ a*^ is independent of proton spin, so that 

(m / |pV<|m/}±,±' = ^(i.i'm/)^/ - l/2|pVa e K |m / - l/2)± + c 2 (, h T'm f )(m f + l/2|pV o£|m, + l/2)±. (29) 
Our results for the required expectation values are: 

/ v A l + 4m2 m 2 Wfjl+Mfl 

(m/l ^ K)± ^' = (2 + ^7+1^ i(7+l)fe+i) J 2 

/ 1 1 + 4m? m 2 f \ 

' ' ± v I tn/o- , ^ **8*J3 ( 30a ) 



(m/|pVo-el TO /> 



(m/lpVCTel™/) 



2 8j(j + 1) j(i + l)(2i + l) 
+2m/ 



/ „ , ^+l)-3m 2 / 

x ^+l)-l/2-m? T ' 2 / +1 ' 

+~2m/ 

± - ± ' ~~ (2£ + 3)(2£+ 1)(2£- 1) 

/ „ ^+l)-3m 2 / 

x \i(i+l)-l/2-m)T' 2j / +1 ' 

(m/|p l pV e |m/)±,±' 



(<M*3+M*0> (30b) 



(^2^-3+^3^-2), (30c) 



(2£ + 3)(2£ + 1)(2£- 1) 

, £(£ + 1) - 1/2 T (* + 1/2) + 3m? \ fe^ + fe<5 i2 



+ 1) T + 1/2) + mj T 



2j + l / 2 



/ , , e(£ + l)-l±(£+ 1/2) -3m 2 f \ 
+ I £{£ + !)- 3/2 ± (£ + 1/2) - m? +' + j £ J fe^s 



(30d) 

Some consequences of (|30|) are 

(m f \p ■ p\m f )±,±, = (m/|m/)±,±' = 1, (31a) 

K-lp-p^Kk.i' = Wkek/>±.±'=±^fy^^y43, (31b) 
(m/lpp-a.Kk.i, = (m/lpp'aJlm,)^^ = 2j(j + 1)(2j + ^ 3 = (2/ + 1)(2j + (31c) 
(to/Ip 1 (p x CT e ) a |m/)±,±< = (m/|p , p i <|m/}±,±'e J - fea = + ^ - e*a3 = + 2 / + i e ' a3 ' ^ 31d ^ 
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where (|31al) and (|31bp serve as checks, and (|31c[) and (I31dp will be needed below. In the reductions of (|31c[) and 
(|31d[) we have applied the useful identity (2j + 1 =p' 1) = + l)/(2/+ 1), and we record here the additional identity 
(2j + 1 =f' 2) = (2j + 3)(2j + l)(2j — l)/(4/(/ + 1)). As always, expectation values of p l p'a^ are obtained from those 
oiffp^a* by multiplication by {p 2 ) n — (m e a/n) . 

The corresponding proton expectation values of p'p'cc are obtained by calculations similar to those used for fffp a^. 
We start from (fT4l) and find 



(m / |p i ^'^|m / >±, ± / = ±'2c(j,±'m f )c(j,T'm f )Re [(m f ~ 1/2$? \m f + l/2)±] , (32a) 
(m/l^o^m/)^ = ±'2c(j,±'m f )c(j,T'm f )Im [(m f - l/2\f p>\m f + l/2)±] , (32b) 
(m/l^oflm/)^ = c 2 (j, i'm/Xm, - l/2|^'|m/ - 1 / 2 )± ~ c 2 (j,t'™/) W + l/2|^*|m/ + l/2>±- (32c) 
The required matrix elements in the orbital plus spin states \n,£,j = £ ± 1/2, raj) of (jlip are 



-m f J(j + 1/2) 2 -m 2 

(m/ - l/2|pV|m/ + l/2)± = V ^ - ^ {(^3 + fa^i) + i (feS,3 + fe^)} , (33a) 



(m / -l/2|pV|m / -l/2)± = 



(j + l/2) 2 + m 2 - m/ \ ^1 + gnjja 
2j(j + 1) ) 2 

' (j + 1/2) 2 - 1/2 - mj + mf 



(m f + l/2\p i p>\m f + l/2}± = 



2j(j + 1) 

' (j + l/2) 2 + m 2 + m/ \ + s i2 6 j2 

2j(j + 1) J 2 

'(j + 1/2) 2 - 1/2 -m 2 - to/ 



6 i3 6 j3 , (33b) 



2j(j + 1) 



&3<Sj3- (33c) 



These lead, through equations (|32|) . to 



. . . T'm f ((j + l/2f-rr 
{mf\p t fa p \m f )±,±> = - - (SnS j3 + (fotfji) , (34a) 

T'm f ((j + 1/2) 2 - to 2 

f \p l p><jl\m f )±,±, = \ ————11 (S i2 5 j3 + S i3 5 j2 ) , (34b) 



™ 2j(j + l)(2j + l) 



(to/|;p*pV||to/)±,± 



±'m/ 



" ' J(j + l)(2j + l) 

(J + 1/2) 2 ~\-m) ±' (j + 1/2)) 6*8. 



((, + 1/2) 2 + mj T ' (j" + 1/2)) ( ^|^ 



Consequences of (|34|) include 



(34c) 



(m f \p-po*\m f )±,±. = {m } W k p \m f )±,±> = ±'^^-S k3 , (35a) 
(m f \p l p-a p \m f ) ± ,±, = + ^ + & 3 , (35b) 

771 f 

(m/|p l (p x (7 p ) a |m / ) ±i± / = T 2 f + i e *° 3, ( 35c ) 

Expression (|35a[) is a check showing consistency with (fTB")) . while (|35b|) and (|35cl) will be needed for the energy shift 
calculations. We note the similarity (and identity when ± = ±') between the electron spin results of (|3 1 1) and the 
proton spin results here. 

The 0{SME ■ a 2 ) energy level shifts coming from SME electron and proton interactions for a state with quantum 
numbers n, £, j — i ± 1/2, / = j ±' 1/2, and nif are 



F 2ep 
& SME 



{{m f W\m f ) ± , ±l (E% + eEfj) + (m f \p^ (a k e Ff jk + o k p eF? k ) |m/) ±1± /} , (36) 
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FIG. 1: Graphs involving photon exchange between the electron (top) and proton (bottom). Part (a) shows Coulomb exchange 
photon that is responsible for the dominant binding. Parts (b) and (c) show SME corrections. In (b) the interaction involves 
the SME coefficients (fci?) Q(i/ 3„ acting between two Coulomb photons. In (c) the SME correction with coefficients (Auf)" acts 
between a Coulomb photon (connected to the proton) and a transverse photon (attached to the electron). 



where e = (m e /m p ) 2 . More explicitly, the 0{SME ■ a 2 ) energy corrections from these sources are: 



F 2ep 
C SME 



5 f£ e | c gM | (2i + 3)(2j-l) I m-f i / n 

+ 27TT ( 27TT ( J 3 + «*) + (^s T> } (2) 2 . (37) 

We find that only four combinations of parameters enter into the final result: c o = tticqq (where cqq = Cu because 
is traceless), c Q = m(cu + c 22 - 2c 33 ), d 3 = mdm + \md 30 - \H\i, and g D3 = -b 3 + m(g W2 - g 2 oi + guo)- (The 
definitions were adapted from Table XVIII of Ref. [8j.) 
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IV. PHOTON-SECTOR ENERGY SHIFTS 



In this section we calculate the energy level corrections at 0{SME ■ a 2 ) coming from SME corrections to the 
propagators for exchange photons. The graphs in question are shown in Fig.[T] The Feynman rules for SME corrections 
to the photon propagator are given, e.g., in the appendix to [22j |. There are two corrections: one (indicated by a large 
filled dot), has the form —2iq a q@ (kF)afj,i3iy, where fi and v are photon indices and momentum q travels from the [i 
to the v index. The second (indicated by a cross) has the form 2{kAF) a £afipv<l^ > where e a(1 ^ is the four-dimensional 
Levi-Civita symbol (€0123 = —1). The counting rules for QED bound states such as hydrogen allow us to estimate 
the orders in a of the graphs shown. Each interaction vertex brings a factor of e oc a 1 / 2 . Three-momenta are of order 
a and energies are a 2 . Wave functions scale as a -3 / 2 in momentum space. The first (Coulomb interaction) diagram 
has an interaction factor (e)(l/<f 2 )(— e) = —4ira/q 2 and a corresponding energy given by 

Ec = I WfUf { ^ m ^' icw) (r|m/>± > ± ' (38) 

where q = p — p' . The scaling laws imply an order 0(a 2 ) for this energy, and indeed the Fourier transform of the 
interaction is ~a/r with expectation value Eq = (—ot/r) = 2E n = —m e (a/n) 2 . We use Coulomb gauge for the 
photon propagator, and take the QED Feynman rules for this non-relativistic system from Non-Relativisitc QED 
(NRQED), whose rules have been tabulated by, for example, Kinoshita and Nio [33| or Labelle (3~3 |. 

The graph of Fig 1(b) has, compared to the Coulomb contribution of Fig. 1(a), two extra powers of the photon 
three- momentum in the SME vertex and an extra Coulomb propagator. It thus has the same order, a 2 , as the 
Coulomb energy, but involves as well the SME parameters {kF)aii(5v- Spatial parts of the Coulomb-gauge prop- 
agator do not contribute at this order in a because they would come with a "dipole" vertex — e(p' +p)/(2m) 
containing an extra power of the three-momenta. The interaction factor arising from the graph of Fig. 1(b) is 
(e)(l/q 2 )i (— 2iq l qi (kF)iOjo) (l/9 2 )( — e ) = —2e 2 {tF)iOjoq l q : '/q' 4 where we have used the antisymmetry of {kF) a ^pu 
in the first and in the second pair of indices. I6| [22l| The extra factor of i before the SME Feynman rule factor is to 
convert to the conventions used by Refs. [33l |34{. The Fourier transform of the momentum-dependent part of the 
interaction term, q % q> /q 4 , is (5y — x l x^ /(8nr), as shown for example in Ref. [35j], so the correction to the potential 
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energy is (a/r)(k F )iojo {x % x 3 — 5ij). This is consistent with the potential found in Refs. [361-138] from the solution of 
the SME-modified Maxwell equations. The corresponding correction to the energy is 



d 3 x{x\m f )* ± ± , (^—(k F )iOjo {x l x 3 -Sij)^j (a|m/)± ( ±/. (39) 



As in our earlier momentum space work, the radial part factorizes, giving (a/r) n = m e {a/n) 2 . We are left with an 
angular expectation value 



m e (k F )iOj(){mf \ [tfx 3 - Sij) |m/)±,±< y—J (40) 



The angular expectation value here is the same in coordinate space as in momentum space, so we can use (I30ap for 
the expectation value of x l x J . We find it convenient to write this expectation value as 

(m f \x^\m f ) ±i± , = (m f \mm f } ± , ± , = ^ - ^/ff' i) ~ (l ~ ftf+Tj) (41) 
where 8fj = SuSji + #j2<^'2 — 25^8^. The resulting expression for the energy shift is 

J 2. (2j + 3)(2j- 1) ( 1 mj \ _ \ faV 
EF = 16,(, + 1) 1^3 - 7(7+1) J ^ j U ' ^ 

where k = m e ((k F ) 1010 + (fc F ) 20 20 + (k F ) 30 3o) and k q = m e ((Mioio + (k F ) 2 020 - 2(fc F ) 3030 ). 

The graph of Fig. 1(c) involves the kA F SME parameters. Because of the Levi-Civita symbol in the 
Feynman rule, only one of the photons can be Coulomb, the other must be transverse. We calculate the 
graph having the transverse photon attached to the electron because it is larger by a factor of m e /m, p 
than the one with the transverse photon connected to the proton. The interaction factor for this graph is 
(-e(p + p') a /{2m e )) (-S ak /q 2 ) i (2{k AF ) i e i0jk q^) (l/q 2 ) (-e) = (ie 2 /m e ) (k AF ye 0tjk q 3 {p + p') k /q 4 . We can write 
(p + p') k = (q + 2p') k — > 2p lk , so the interaction factor becomes (— 2ie 2 /m e ) (k AF yeij k q J p' k /q A . This interaction 
factor is not just a function of the relative momentum q, so we work out the Fourier analysis in more detail. The 
energy shift from Fig. 1(c) is 



2ie 2 /, m [ d3 P d3 P' fl j P' k 



Eaf = — {kAFy ^ k J ( 27r )3 (2ff) a <?K>3:,±' [-^r) tf'\ m f 

2'' " . , [ d 3 p '/•'/'' ,:; ,:; ■ ( <l' 



;±,±' 



m e 



(k AF y^ k J ^0^d 3 xd 3 y(x\m f )^ ± , e#* (*\ *V*e-^%|m / ) ± ,±* 



m e 
-2ie 2 



(k AF fe ijk J d3xd3 y^30j3 (2\m f )i, ± , e -#'-V-*> (-iVj) (y \m f ) 

(k AF ye ijk J<Px{S\m f y ±i± , (-<v£) (2|m/)±,±'- (43) 

So this energy is just proportional to the expectation value of the orbital angular momentum L = x x p. We use the 
Wigner-Ekhart theorem as in (|16p to write L — > inside the |m/)± ±/ expectation values where 

^ (21 +l)(2j + 1) 2 (2/ + 1)' 1 ' 

We found ^ by direct evaluation using the states |to/)± i ±' of (fM]) . and checked using £^ + £ e /2 + £ p /2 = 1, which 
follows from L + a e /2 + a p /2 = F. The E AF energy shift is thus 

E AF = ^- (~) (m f \(k AF yL l \m f ) ± , ± , = (k AF ) 3 ^m f . (45) 
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V. RESULT, AND APPLICATION TO THE 25 IS TRANSITION IN HYDROGEN 

The total energy shift at order 0(SME ■ a 2 ) is the sum of the fermion-sector shift E s e ^ E of (|37|) and the two 
photon-sector contributions Ep of f4"2")l and Ej^p of f4"5)) . It is 

Elu* = { - I (clo + «*,) - f ft + W + ^ - 1) A _ , ^ + £?q _ ~ KQ 




2/ 

The SME coefficients appearing here are those for the lab frame. The SME coefficients in the lab are related to those 
in the conventional sun-centered inertial frame by a Lorentz transformations involving both a boost and a rotation. 
The expression for E 2 ME given by (|46p is our main result. 

The transition between the 25 and 15* states is of particular interest because it has an anomalously small intrinsic 
broadness and it has been measured with great precision. The 25 — > 15 electric dipole decay into a single photon is 
forbidden by the usual Ai = ±1 selection rule. The decay happens b y a slower two-photon transition with a lifetime 
of about 0.12s (as opposed to the 1.6ns lifetime of the 2P state). 1391 . Eol] The corresponding uncertainty-relation 
width of the 25* state is about one Hz. A recent measurement gave [2l[ 

hs-2S = 2 466 061413187035(10)iJz (47) 

for the transition frequency with a precision corresponding to a fractional uncertainty of 4.2 x 10 -15 . 

We can easily use our results to evaluate the SME corrections to S state energies and splittings. For 5 states the 
momentum product p l p> is effectively ^5ij and the expectation values of (f30|) and (j34|) become 

(m f \p'p j \mf) + ,±> = ^Sij, (48a) 

(mf\fp>a*\m f ) + ,±> = ^ij^fem/ =^8 ij 8 k3 m { , (48b) 

{m f \p l p J (jp\m J ) + ±, = ^5ij^°S ks m } = ^5 i:j 5 kz mf (48c) 

where ± — >■ +1, j — >• 1/2, and ±' — > +1 for states with / = 1 (which states are the only ones with non- vanishing m/). 
It follows from (JTUJ) , Q, (@0j), |45]), and (|48]), or directly from (J46J), that the SME energy correction at 0(SME- a 2 ) 
is 



R 12 

fi SME 



{~l (coo + <o) - \*o + ^ (di + e% - gU ~ ^3) } (49) 



for a state with a particular value of n. The 25 — 15 transition in hydrogen is dominated by two-photon emission 
that is subject to the selection rules: Af = 0, Am/ = 0. [4l[ It follows that there is no SME contribution to the 
25 — 15 energy splitting at 0(SME ■ a ). The leading correction to the splitting is at 0(SME ■ a 2 ): 

AE SME (2S - IS) = j ~ (c% + ec p 00 ) + 2k Q - m f (d e 3 + e% - g'm ~ ^3) } (50) 

These results are consistent with Eg ME —> — (| c oo) °f Altschul [IH, who discarded coefficients other than c e 

in J2|) as already well-bounded and did not consider the k F contribution, and with Av — > — for the difference 
between the / = l,m/ = 1 and the / = 0, rri f = frequency shifts of the 25 — > 15 transition, calculated as an 
example by Bluhm, Kostelecky, and Russell (12f . 



VI. DISCUSSION 



Having values for the 0{SME ■ a 2 ) corrections to the hydrogen energy levels makes it possible to consider tests in 
hydrogen for additional SME parameters besides just B k and B v k . In fact, B k and B v k cancel entirely for the 25 — 15 
energy difference, leaving the 0(SME ■ a 2 ) contributions as the leading SME corrections. In addition, it would be 
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easy to find the corresponding corrections for anti-hydrogen, since the SME coefficients for antiparticles are related 
in a simple way to the coefficients for the corresponding particles. Specifically, the coefficients for positrons and anti- 
protons are the same in magnitude as the coefficients for electrons and protons, but the a M , d^ v , and coefficients 
change sign. [13 . [l8| Comparisons between hydrogen and anti-hydrogen have the potential to probe separately the 
various parts of Bk and to more easily distinguish between the effects of dk, which changes sign for anti-hydrogen, 
and coo, cq, an d <)Dk, which do not. The effects of the proton parameters are suppressed in hydrogen by a factor 
of e = (rn e /m p ) 2 w 3 x 10~ 7 relative to electron effects because the typical proton speed squared in hydrogen (v 2 ~ 
(\p\/m p ) 2 ~ (am e /m p ) 2 ~ 1.6 x 10~ n ) is so much less than that of the electron (v 2 ~ (\p\/m e ) 2 ~ a 2 ~ 5 x 10~ 5 ). 
Much larger proton (and neutron) effects can be found in systems with a nucleus containing more than just a proton, 
as nucleon speeds in nuclei are much larger, typically v 2 ~ 10~ 2 . fl9| 

We expect that there are additional energy corrections of 0(SME ■ a 2 ) in hydrogen due to interactions such as 
an SME fermion line correction of the form A + B ■ a combined with exchange of, say, a transverse photon. Such 
un-calculated corrections depend on the same SME parameters as in (|17p but are higher order in a. We have restricted 
our attention to the leading order effects of any particular SME parameter. 

In our work up until now we have assumed the existence of a small external magnetic field. This field has the 
effect of breaking the degeneracy over mt so that the states |m,/)±,±' of (TH)) are non-degenerate and non-degenerate 
perturbation theory can be used in our calculations. We note that the direction of this external field determines which 
component of the relevant SME parameters can be probed by the hydrogen experiments. The situation changes if this 
field is completely absent or shielded so well that its effect is less than that of the SME. In that case, the degeneracy 
over m/ would be broken by the 0(SME ■ a ) perturbations, leading to a lowest order SME energy shift of 

£ S °me (m f \ {A e + AP + B%al + B p k a k p ) \m f ) ± # = (A e + A p ) + \^ e B e + t p &\ m f , (51) 

where now the z axis is determined by the direction of £, e B e + (,pB p '■ The 0(SME ■ a 2 ) corrections take the same 
form as (|46|) except that now the z axis is determined by the SME vector £, e B e + £, p B p instead of an external magnetic 
field. 

We note that frequency measurements are always done relative to some standard. For example, in the 2S — > IS 
measurement of Parthey et al. [2l| , the measurement was done relative to a cesium fountain atomic clock. A SME 
effect on hydrogen is in principle only detectable if it affects hydrogen and the reference clock differently. In our case, 
the main SME effects for hydrogen are through B k , B k , the photon coefficients ho, kq, and (kAF) k , arid the electron 
coefficients c^q, Cq, d%, and g e Dk (assuming that the tiny e factor causes the proton contribution at 0{SME ■ a 2 ) to 
be negligible), while nuclear effects are important in cesium. "Clock comparison" experiments have been discussed 
extensively by Bluhm et al. in the context of the SME. [29|, |30| 
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Appendix A: Effect of a coordinate transformation 



Many SME coefficients or combinations of coefficients are unphysical and cannot be observed. Unobservable fermion 
coefficients, such as a M + me M , Cp V — c VIJ ,, and f p , along with the corresponding field redefinitions used in the demonstra- 
tions of unobservability, have been discussed by many authors. 0, [25l - l27j Invariance of the physics under particular 
coordinate transformations can be used to obtain relationships among certain photon and fermion coefficients, specif- 
ically between a subset of the (kF) pupi7 photon coefficients and the symmetric and traceless c M „ fermion coefficients, 
[tl l36l I371. |42M44| Suppose we consider a theory containing only electrons and photons, and only the c M „ and {kp) lliV p a 
SME coefficients: the Lagrangian reduces to 

£ = W + c» v )D v il>- - -F^ {rjTrr + (k F m F pa . (Al) 
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The nineteen inde pen dent coefficients of (k,F) puf " T have the symmetries of the Riemann tensor and are double traceless 
by convention: [(I |36| 



(&f) pupa — V*F) papv — {kF^vppa \ 
+ (k F )pavp = 0, 

(fc F rV = 0. 



(A2a) 
(A2b) 
(A2c) 



The {kF)^ upa tensor can be decomposed into a ten-component completely traceless part C^ vpa and a nine-component 
part defined in terms of a traceless tensor k^ = (kF) pau ' a according to: (38l. lisl l46j 



We restrict our model further to have C pvpcr = 0. The coordinate transformation 



1 - 



x>- -> x 



lr, 



9a? ' 



(A3) 



(A4) 



has unit Jacobian due to the tracelessness of k p ,, and allows us to write C as 



lr, 



C = ^% tT + - -hT pW - - - A F'rf p rf°F' 



pa 



(A5) 



where ip' = ip(x'), D' v — d' v + iqA' u (x') : etc. So a theory of electrons with the usual SME c^ v coefficient and 
photons with the SME (kp)pvpa coefficient given by the k^ v part of (|A3[) is equivalent to a theory of electrons with 
Cfj, v — > Cfu, — (l/2)kp I/ and a conventional quadratic photon Lagrangian. This equivalence must be evident in the 
electron-photon part of our result for hydrogen energy levels. However, we do not expect the energy in the unprimed 
coordinate system to equal the energy in the primed coordinate system-instead, it should be related according to 
transformation (|A4[) . For a state with energy E and time dependence exp(— iEt) in the unprimed system, the time 

dependence in the primed system takes the form exp iE(l + koo/2)t'^ . The corresponding primed system energy 
will be 



E' = ( 1 + ^fcoo 



E. 



(A6) 



The relevant part of our hydrogen energy result (1461) is 

2 



E = m jl - c q - 



15 2 



15 2~ ca 
77 - 77C00 + 77^00 + Xdl 



ij I °ij 



Ik 



(A7) 



where A is a relatively complicated function of the angular quantum numbers and (|A3|) with C p,I/f " 7 = leads to 



(kF)iojo = 77 [hj - Sijk, 



Urn 



(A8) 



so that Sij(kF) i0 jo = (kF)aOaO = —koo and Slj(kF)iojo — ^S^kij. We transform to the primed system by first making 
the replacement k^ — > followed by the replacement — > c^ v — (\/2)h^ v in (|A7|) . and obtain 



77 - 7: ( coo - 77^00 ) + X8% ( c 



h 

2 13 



(A9) 



To first order in SME coefficients this has the required form (|A6 
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